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^\ ' We discuss a way to deduce parameters of accreting black holes. The method employs properties of spectral 

^^ ■ features observed in radiation from an accretion disk. It is applicable to sources which exhibit periodic modulation 
of variability, provided: (i) Gravitational field is determined by the black hole and described by the Kerr metric; 
(ii) A thin accretion disk of a negligible mass lies in the equatorial plane of the hole; (iii) A secondary object (also 
with negligible mass) moves on a slightly inclined almost circular orbit around the black hole and passes periodically 
through the disk; (iv) The collisions result in observable photometric and spectroscopic features (temporal variability 
of the radiation flux and of spectral- line profiles produced in the disk) , which show frequencies of orbital motion 
and of latitudinal oscillations; (v) One can measure width of the spectral line from hot spots arising in the disk 
due to collisions with the orbiter, and/or detect the predicted low-frequency oscillations, induced in the disk. 
These frequencies and the line width provide enough information to determine in physical units three parameters 
characterizing the source — the mass and angular momentum of the central black hole, and radius of the orbit of 
3 . the secondary. 
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1. Introduction 

There is an ever-growing theoretical and observa- 
tional evidence that the dark masses present in some 
galactic nuclei and X-ray binaries are black holes 
(Rees 1998). Standard interpretation of observational 
data employs a black hole and an accretion disk or a 
torus (Kato et al. 1998). A number of ways have been 
suggested of how to deduce the properties of putative 
holes. The first rough estimates of their mass were 
based on energy considerations and limits implied by 
the shortest time-scales of variability. More precise 
methods became possible with modern observational 
techniques like HST, VLBI and X-ray satellites. Our 
present paper is relevant for sources with a rotating 
black hole surrounded by a thin accretion disk in the 
equatorial plane and a much less massive secondary 
orbiter passing periodically through the disk. First 
we briefly recall relevant spectral properties which 
have been discussed as possible observational signa- 
tures of a black hole (we do not mention other ap- 
proaches which involve different physical mechanisms, 
e.g. gravitational waves). 

In galactic nuclei, the existence and the size of 
the dark central mass are deduced from the profile 
of the surface brightness of the nuclei and from the 
spatial distribution and dynamics of surrounding gas 
and stars (Kormendy & Richstonc 1995). Current 
optical and radio studies still probe only scales above 
10 4 gravitational radii of the nucleus which is insuf- 
ficient to resolve imprints of relativistic effects in the 
innermost regions around the centre, in particular to 
deduce the rate of its rotation. However, X-ray ob- 
servations of the fluorescence iron Ka line (Nandra 
et al. 1997), especially its broad and variable profile 
skewed to lower energies, indicate that the Dopplcr 
and gravitational redshifts play a role and the em- 
mision takes place very near a rotating black hole. 
Several ways have been proposed to estimate param- 
eters of the black holes in active galactic nuclei which 
rely on the source variability: for example, limits on 
the black-hole mass and its rate of rotation can be 
imposed from time delays between variations of the 
emission-line strength and of the continuum, and from 
temporal changes of the observed emission lines (e.g., 
Blandford & McKee 1982; Stella 1990; Bromley et al. 
1997). 

The evidence for stellar-mass black holes comes 
from observations of Galactic X-ray binaries where 
they are indicated by low luminosity of one of the 



components, temporal variability, spectral features 
(mainly the presence of relatively strong ultrasoft 
component and of the hard X-ray tail), and the mini- 
mum mass of the dark component, as estimated from 
the mass function (e.g. Charles 1997). It has also been 
argued (Narayan et al. 1997) that the black-hole X- 
ray binaries could be distinguished by larger variation 
in luminosity between the bright and the faint states 
than is expected in the sources with neutron stars. 
This way of identifying black holes follows from the 
low radiative efficiency of advection-dominated accre- 
tion flows around black holes. Present observational 
evidence offers several objects of this type which show 
relatively stable periodic modulation (van der Klis 
1997). 

In the present contribution it is assumed that the 
disk remains in the equatorial plane and the mod- 
ulating source is represented by a blob which keeps 
its identity for several orbital periods. A specific 
scheme for the signal modulation can be described 
in words as follows: The orbiter intersects the disk 
periodically and pulls gaseous material out its plane. 
This material temporarily obscures the disk, mainly 
at the radius of intersection, and affects radiation flux 
and emission-line profiles formed in that region and 
characterized by the observed widths. Phcnomeno- 
logically, a spot forms in the disk and then it orbits 
around the center, surviving several orbital periods. 
Radiation of the spot modulates the observed radi- 
ation at the Keplerian frequency which is different 
from the frequency of collisions with the disk when 
the central black hole rotates. Both ingredients of 
the model, i.e. local physics of collisions with the disk 
and observed signal from orbiting spots, have been 
described in literature (recently, Ivanov et al. 1998; 
Karas 1997). In usual terminology the blob repre- 
sents a hot spot, however, any perturbation of the 
disk surface emissivity (either bright, or less promi- 
nent due to obscuration) suffices for the purpose of 
phenomcnological description. 

Subsequent paragraphs describe technical details 
of the above-mentioned scenario. We will also show a 
predicted variable spectral feature involving relevant 
time-scales in a simplified situation. 

2. Assumptions 

Consider a system of a rotating black hole, an equa- 
torial thin disk and a secondary orbiter (which may 
be a low-mass black hole -- see Syer et al. 1991; 



Vokrouhlicky & Karas 1993). If both the disk and 
the secondary have negligible masses with respect to 
the central hole, the gravitational field is deteremined 
solely by the central black hole. In this situation, 
motion of the secondary is governed by Kerr metric, 
which in usual notation (Misner et al. 1973) acquires 
the form 



ds 2 = -AXL/TMr + HA-Mr^ + Edfl 

+AT 1 - 1 sin 2 e {d(j> - u K dtf ■ 



(1) 



here Boyer-Lindquist spheroidal coordinates (t, r, 8, (j>) 
and geometrized units (in which c = G = 1, c be- 
ing the speed of light in vacuum and G the gravi- 
tational constant) have been used, M and a denote 
mass and specific rotational angular momentum of 
the centre, A = r 2 - 2Mr + a 2 , £ = r 2 + a 2 cos 2 6», 
A= (r 2 + a 2 ) 2 - Aa 2 sin 2 6, and uj k = 2Mar/A. 

World-line of the secondary is very close to a 
geodesic in the Kerr field, provided that the dynam- 
ical effect of passages through the disk is only weak, 
that tidal interaction is negligible (the secondary is as- 
sumed much smaller than the typical curvature radius 
of the field around), and that gravitational radiation 
is ignored. The trajectory undergoes three secular 
changes due to the interaction with the disk (Syer 
et al. 1991; Vokrouhlicky & Karas 1993): long-term 
decrease of the semi-major axis (spiralling towards 
the centre due to energy dissipation in collisions with 
the disk), circularization (the orbit becomes spherical, 
r = const), and gradual tilting of the orbit into the 
equatorial plane (the orbit declines to the disk). Since 
the time-scale for circularization is shorter than, or of 
the same order as, the time-scale necessary to drag 
the orbit into the disk, one can assume that the sec- 
ondary follows a nearly equatorial spherical geodesic 
at late stages of evolution of the hole-disk-secondary 
system. The secondary remains on this type of orbit 
for a relatively long time. Let us note that we ignore 
gravity of the disk itself; the case of a massive disk 
was treated recently by Vokrouhlicky & Karas (1998). 

The resulting system is characterized by two angu- 
lar frequencies — that of the azimuthal revolution of 
the orbiter around a Kerr black hole, u>, and that of 
its latitudinal oscillations about the equatorial plane, 
iloo. Both frequencies are in principle measurable at 
infinity provided that the passages of the orbiter pro- 
duce strong enough modulation of the disk radiation 
(Karas & Vokrouhlicky 1994 illustrated, by Fourier 
analysis of simulated photometric data, how these two 
peaks can be recognized in the power spectrum). 



Relevant parameters of a quasi-circular orbit pro- 
cessing with low amplitude about the equatorial plane 
of the Kerr spacetime are given in next section. Rela- 
tions derived in Sect. 4. provide the parameters of the 
central black hole in terms of four observables — the 
three frequencies (f2oo, w+, and k, as they are intro- 
duced in the text) and a spectral-line width. Equa- 
tions can be solved easily if the secondary is not too 
close to the black hole (cf. the Appendix). 

3. Parameters of the Precessing Orbit 

For a general bound geodesic in the Kerr space- 
time, the frequencies of the azimuthal and latitudi- 
nal motion can be given in terms of elliptic integrals 
(Karas & Vokrouhlicky 1994). Relevant expressions 
acquire a simpler form in case of spherical geodesies 
(Wilkins 1972) and simplify still further for almost 
equatorial orbits. The azimuthal angular velocity 
lu = dcf>/dt (with respect to a distant observer at rest) 
is then approximated by the Keplerian circular fre- 
quency 



uj± = (a+ l/y±) 1 , 



(2) 



where y± = y(u>±) = ±\J ' M /r 3 and the upper/lower 
sign corresponds to the prograde/rctrograde orbit. 
We keep both cases for completeness, but only pro- 
grade trajectories (plus sign) will be considered later 
(the accretion disk is more likely to be corotating with 
the central black hole and the interaction with the 
disk makes also the secondary corotate eventually). 

The proper angular frequency |f2| of small lat- 
itudinal harmonic oscillations about the equatorial 
plane is given, for a spherical orbit with steady radial 
component of acceleration, by (de Felice & Usseglio- 
Tomasset 1996; Scmcrak & de Felice 1997) 

n 2 = (u'/r) 2 {Au 2 + 2y 2 ± [a - (r 2 + a 2 ) u] 2 } . (3) 

Here u> = const, and the time component of four- 
velocity, given by (u*)~ 2 = -g tt - 2g t(j> u! - g<j><j>u 2 , is 
approximated by the equatorial value 



(u 1 )- 2 = 1 - (2M/r)(l - aujf - (r 2 + a 2 ) < 



(4) 



Frequency of the collisions with the disk is twice the 
proper frequency |fi|; the corresponding value mea- 
sured by a distant observer is |fioo| = |f2|/u*. From 
eqs. (2)-(3), the explicit expression for a, free orbit is 



<>: 



w | 



(i 



4at/± + 3a 2 /r 2 ) 

+ 2(VAfrTfl) 2 



(5) 



We remind that the frequency of latitudinal oscilla- 
tions stands in the dispersion relation for corrugation 
waves in accretion disks. It has been suggested that 
fioo could be detected in some QPO sources (Kato 
1990). 

The two eqs. (2) and (5) are not enough to fix the 
three unknowns M, a and r. Another piece of in- 
formation can be provided by the width of the emis- 
sion line which is formed and modulated at the ra- 
dius of successive collisions between the secondary 
and the disk. A stationary disk produces the well- 
known double-horn line profile which is more or less 
pronounced according to the inclination angle of the 
source. Radiation from each element in the disk ex- 
periences a frequency shift (Fanton et al. 1997) 
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where ei is an emission direction cosinus (an az- 
imuthal component of the unit vector along the di- 
rection of emission of a given photon, measured in 
the emitter's local frame), and 
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(7) 



follows from (4) and (2). The observed total width 
of the line arises from the different frequency shifts g 
carried by photons which reach the observer at infin- 
ity. Since g depends on the emission cosinus e;, the 
maximum width of the line Sg (which is a result of 
integration over one entire circle) corresponds to the 
range of e; consistent with the escape to infinity, Sei. 
We have from cq. (6) 



(8) 



(6g) = (6ei)u t u ± VA 



and thus 



(Sg) 2 _ M l-2M/r + a 2 /r 2 
(Se^) 2 ~ V 1 - 3M/r + 2ay± 

M A 

r A- (VMr=Fa) 2 ' 



(9) 



From this equation it appears convenient to include 
the parameter 5 in our consideration because it can 
be written in the form independent of inclination an- 
gle. To determine 5 in terms of Sg, one needs to fix 



the range of Sei. This can be easily estimated for 
r ^> M (cf. Appendix). In general, Sei is at most 2 
(which is acquired for 8 = ir/2, i.e. edge-on view of 
the disk), and it decreases with the inclination angle 
O of the black hole-disk system relative to the line of 
sight (it also depends, quite weakly, on the rotational 
parameter a and the radius of emission r) . It appears 
unpractical to deal with a lengthy analytic expression 
for Sei. Instead, given a/M, r/M and O , the actual 
value of Sei can be found by numerical ray tracing 
(Fanton et al. 1997). 

We simulated a time-variable line profile to illus- 
trate the main features expected in the observed sig- 
nal (we used a modification of the code described 
by Karas and Vokrouhlicky 1994). Figure 1 shows 
the measured count rate in the line (background sub- 
tracted) as a function of energy and orbital phase in 
the equatorial plane of an extreme (a = M) Kerr 
black hole. In this example, only that contribution 
to the total light is taken into account which orig- 
inates at 10M < r < 14M where an orbiter affects 
the disk. Energy is normalized to the emission energy 
of the line in the local frame of the disk material, as 
usual (energy axis thus indicates the redshift factor 
g). Phase axis is normalized with respect to Keple- 
rian orbital motion at the central radius r = 12 M, 
and interval of two periods is shown for clarity. One 
can observe: 

(i) Modulation of the count rate in the line by the 
orbiter crossing the disk upwards at radius r with 
period T = 27r/£!oo(r). A flare is triggered at the 
moment of crossing and then falls off with time as 
the orbiting spot decays and the material swept out 
of the disk cools down (characteristic time-scale for 
this decay is taken to be two orbital periods; it is a 
phcnomenological parameter here). 

(ii) Another modulation by the orbiting spot with 
period t = 2tt/lo + . Corresponding maxima of the 
count rate can be seen at phase equal to 0.5, 1.5, . . . 
(due to Doppler boosting). 

(hi) Orbital variation Sg of the redshift factor (also 
with periodicity t) in projection down into the energy- 
phase plane. 

In Fig. 1 inclination angle has been taken 40 deg. 
Obviously, any non-zero inclination would lead to a 
qualitatively similar picture. However, at higher in- 
clinations (> 60 deg) another local maximum would 
develop in the observed flux due to the lensing effect. 
(Its position is advancing in phase by w 0.25 with 



respect to the Doppler peak.) Difference in the two 
frequencies, ftoo and w+, discussed in previous para- 
graphs translates to the time gap T — t. Let us em- 
phasize that this figure offers a qualitative illustration 
of the predicted time-variable spectral feature. In a 
more realistic model one needs to take into account 
a proper method of background continuum subtrac- 
tion (underlying continuum has been ignored here). 
Mutual interplay between spectral features (around 
6.4 keV) and continuum has been subject of recent 
works (Young ct al. 1998; Zycki et al. 1998; Martoc- 
chia et al. 1999) which show substantial smoothing of 
the profile in some situations. 

Uncertainties connected with the badly known 5e ? 
can be overcome by introducing another independent 
observable. A promising possibility arises from a more 
realistic view of the star-disk collision. As mentioned 
above, the collisions disturb the disk periodically with 
latitudinal frequency (5). If the orbiter passes across 
the disk at small r, it may induce waves in the disk 
matter which remain trapped in the range of several 
gravitational radii and are governed by epicyclic fre- 
quency k (Okazaki et al. 1987; Nowak & Wagoner 
1992): 



u4 (1 - 6M/r + 8ay± - 3a 2 /r 2 ) 



A 



(10) 



These modes can be excited and trapped in a rather 
narrow range of radius, and one can thus expect that 
they will reflect just the orbital radius of the satel- 
lite. However, the problem of distinguishing between 
different possible modes is a subject of much debate 
and has not yet been settled. 



4. Determining M, a, and r 

The formulas (2), (5), (9) and/or (10) provide a 
closed system of ordinary equations for M, a and r 
(and possibly 5ei, if one succeeds in measuring all 
four observables) as functions of the quantities u)±, 
|fioo|, \8\ and/or |k|. The equations can be tackled in 
different manners depending on which of the quanti- 
ties is determined with the best confidence. Now we 
discuss derivation of the most desired combinations 
of parameters, r/M and a/M, from different starting 
relations: 

(i) With the knowledge of |ftoo| and |k| one can start 



from cqs. (5) and (10) and find 
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which can be solved with respect to M/r and a/r 
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(13) 
(14) 



The values of r/M and a/M are reached then by ob- 
vious manipulations. One more relation is needed to 
obtain the absolute values M, a and r. Using eq. (2), 
we fix 



r \r V M 



w±; 



(15) 



M and a follow from above. 

(ii) Starting from eq. (9), a convenient possibility is 
to combine it with cqs. (5) and (10). We find 



S 2 = 



2M2ft? 



3r ft? 



(16) 



and thus the value of M/r follows immediately. This 
can be compared with the result (13) to fix (or verify) 
the value of 6ei. 

(iii) Suppose we know u>±, |ftoo| and \5\. Introducing 
a = u ± 1 -y ± 1 (17) 

from eq. (2) to eqs. (5) and (9), we obtain 

2 _ 3(l-^±/y ± ) 2 



4w±y± + ft 2 



5wl 



(18) 



y± is then given by quartic equation. This is rather 
cumbersome to be handled analytically and it is more 
suitable to find the solution numerically for each given 
set of data. However, a simple explicit solution can 
be found analytically if r 2 ^$>a 2 ; see the Appendix. 

(iv) The case when one measures ui±, \k\ and |<5|, 
and thus turns to eqs. (2), (10) and (9), also leads to 
a quartic equation. 
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Fig. 1. Simulation of the time-variable line profile from a hot spot orbiting with the disk. A flare occurs at 
the moment of the star-disk collision (the sharp increase in flux), giving rise to a spot on the disc surface. 
Observed count rate is then modified by orbital motion of the spot (period t) and by relativistic lensing (with 
a corresponding change of energy in the range Sg). The spot decays gently with time, until another one is 
generated in a subsequent collision (i.e. after period T). This graph has been obtained by ray tracing in the 
Kerr metric. Decay time of the spot is taken as a free parameter. See the text for details. 



In the Schwarzschild case, a = 0, eqs. (2) and (5) 
are reduced to 



u± 



= 2 /±=±|n oo |=±VMA 



eq. (9) reads 

X 2 2 2 * - 2r M 

o = r uj i — -77- 

± 1 - 3r 2 c4 
and eq. (10) 

k 2 = w±(l-6M/r). 



(19) 



(20) 



(21) 



The physical solution of eq. (20) is 

~2 



(4w 



I)- 1 



35 2 + 1 - v/(3<5 2 + l) 2 - 8(5 2 (22) 



while that of eq. (21) 

r 2 = (6C4)- 1 (1 - « 2 /w±) 



(23) 



M follows then from eq. (19). 

Note that the above formulas are written in ge- 
ometrized units. Corresponding quantities in physical 



units are obtained by the following conversions: 
MP h ? s _ M 

M* s ~ 1.477 x 10 5 cm ' 



a phys = ca, 



r phys _ 



also 



a o 



phys 



phys 



M GMp^/c' M GMp^Vc 2 



(24) 
(25) 

(26) 



To obtain frequency in physical units [Hz] (either from 
u> or f2), one uses the relation 

,.,phys rllJ 

./•phys = = — = (4.77 x 10 9 cm/s)u>. (27) 

2tt 2tt 

In usual notation of relativistic astrophysics the (ge- 
omctrized) frequencies are scaled by M _1 (as in 
Tab. 1). Their numerical values must be multiplied 
by factor 

M = (3-231 x 10*) (— j [Hz]. (28) 

5. Discussion 

We examined an interesting possibility of deter- 
mining the intrinsic parameters of a black hole - ac- 
cretion disk system in sources with periodic modula- 
tion of variability which is caused by an orbiting satel- 
lite. The method combines different pieces of infor- 
mation contained in time- variable spectral and photo- 
metric signal, which have been discussed (widely but 
separately) in the recent literature. It can provide all 
the relevant quantities in physical units, but of course 
it can also be combined with another independent de- 
termination of some of the parameters. We assumed 
that the orbiting secondary is only weakly affected by 
the disk, which allowed us to approximate its orbit 
by a nearly equatorial spherical geodesic. This re- 
striction must be abandoned if the secondary is not 
compact enough. 

The proposed approach can be applied to galactic 
nuclei as well as to stellar-mass objects in the Galaxy, 
but identity and physical properties of the secondary 
depend on the mass of the centre around which it 
revolves. In the latter case the lifetime of the sys- 
tem is more restricted by tidal interactions and the 
secondary appears unlikely to survive long enough. 
Effects of tidal distortion of a stellar body passing 
near a black hole have been discussed by several au- 
thors who show that the star becomes squeezed and 



heated (Rees 1988; Luminet & Marck 1985; Evans & 
Kochanek 1989). Tidal disruption occurs in case of 
a close encounter, when the satellite plunges below 
the tidal radius R t ~ 10 11 (M/M*) 1/3 {R*/R Q ) cm, 
where M* and i?„ denote the mass and the radius 
of the satellite. On the other hand, distortions of a 
compact satellite near a supermassive black hole are 
negligible. Duration of a spot or a vortex in an accre- 
tion disk has been also discussed by several authors 
but remains still uncertain (Adams & Watkins 1995; 
Bracco ct al. 1998). We assumed that the lifetime 
exceeds the corresponding orbital period. 

Several possible targets have been discussed al- 
ready. For example, optical outbursts in the blazar 
OJ 287 have recently been modelled in terms of a 
black-hole binary system by Villata et al. (1998). The 
source exhibits several time-scales: feature-less short- 
term variability, 12-yr cycle, and, possibly, a 60-yr 
cycle. These authors, however, presume both compo- 
nents to be of comparable masses while in our calcula- 
tion frequencies are determined under the assumption 
that the secondary is much less massive than the pri- 
mary (cf. also Sundelius et al. 1997). There are more 
BL Lac objects which show a periodic component of 
optical variability but time-scales are always > lOyrs 
in the optical band, and statistical compilation of data 
is therefore very difficult (Liu et al. 1997). As a better 
example, a 16-hr periodicity was reported in the X- 
ray signal from the Seyfert galaxy IRAS 18325-5926 
(Iwasawa et al. 1998). 

We conclude with M. J. Rees (1998): There is a 
real chance that someday observers will find evidence 
that an AGN is being modulated by an orbiting star, 
which would act as a test particle whose orbital pre- 
cession would probe the metric in the domain where 
the distinctive features of the Kerr geometry should 
show up clearly. Indeed, to be able to resolve the ef- 
fects of general relativity in AGN, the satellite would 
have to orbit near the centre and produce periodicity 
on a few-hours time-scale. 
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Table 1. Accuracy of approximative relations (30)-(31). 
Here, a and r stand as input parameters of the sys- 
tem (black-hole angular momentum and radius of the 
orbit of the secondary, respectively). Corresponding 
values of a and f are deduced from our approximative 
equations as given in the Appendix. 
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0.08711 


0.08279 
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0.03065 
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40 
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0.1596 
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40.3 
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6. Appendix. A Simple Explicit Solution of 
Eqs. (2), (5) and (9) for r 2 ^>a 2 

It is relevant to suppose that the secondary orbits 
at r > 10M, and typically at r <~ 30M. In this region 
a significant part of the disk radiation arises. Also, 
the secondary would be tidally disrupted if it plunges 
too close to the black hole horizon (to r ~ M), while 
at too large orbital radius the relativistic dragging 
effect vanishes (with ex r~ 3 ) and the frequencies u>±, 
|fioo| would become indistinguishable. 

At large enough radii (r ~ 30M), several terms 
in eqs. (2), (5) and (9) become negligible when com- 
pared with the rest. Ignoring them, one arrives at a 
simplified set of relations which keeps an acceptable 
precision of the resulting parameters without the need 
to solve the fourth-order equation. (Exact equations 
can be used to improve the estimates numerically.) 
One chooses some particular approximation accord- 
ing to what precision is desired at each of the quan- 
tities M , a, r, and also according to precision of each 
of the input data |u>±|, |f2oo|, \S\. For instance, since 



a V± < 1/30, we can take 

w±-y±(l-aj/±)> y± ~ w±(l + ow±) 

(for the determination of M from the known r, or vice 
versa, it is even sufficient to neglect ay± altogether 
and use just the Schwarzschild form of equations). 
Other two relevant bounds are: 



a 2 /r 2 1_ 
2M/r < 20' 



a 2 /r 2 

ay± 



gy± 

M/r 



An example of a reasonable estimate is provided as 
follows (we suppose a prograde orbit of the secondary, 
and consequently we take plus sign in the relations): 
Restricting to the Schwarzschild limit of eq. (9), 



S 2 = 



M 1 - 2M/r 
~V l-3M/r' 



(29) 



one comes to an approximative value for M/r (the 
approximative values will be denoted by a tilde be- 
low), 



AM/f = 1 + 3S 2 - a/(1 + 3<5 2 ) 2 - 8S 2 



(30) 



Approximative expression for a/r follows then from 
eq. (5): 




-r • ( 31 ) 



Finally, eq. (2) yields 1/f according to the relation 
(15). The unknowns M, a/M and f/M are reached 
by obvious combinations. 

The results of this approximation are illustrated in 
Table 1 for several typical a and r, together with the 
respective values of the observable quantities. Here, a 
and r are assumed to be given (in units of M; see the 
text above for conversions to physical units). For a 
prograde orbit of the secondary, lu + and ri oc (oj + ) are 
the corresponding azimuthal and latitudinal angular 
frequencies, 5 = 5(lo + ) is the dimcnsionless param- 
eter defined by eq. (9). From these quantities one 
computes the estimates a and f (in units of M) using 
our approximate eqs. (30)-(31). Notice that accept- 
able results are reached for r > 10M where approx- 
imative values are close to the simulated exact ones. 
Accuracy decreases if the secondary orbits too close. 
On the other hand, too distant orbits are also unsuit- 
able: dragging effects weaken quickly with distance 



from the centre and the values of uj + and 1]^, be- 
come too close to each other. In other words, the 
minimum sampling rate (inverse of the Nyquist crit- 
ical frequency; Press et al. 1992) is restricted by the 
difference r = f^ 1 — ui+ , requiring that both fre- 
quencies can be safely distinguished in the observed 
signal. For example, for M = 10 8 Mq, maximum ro- 
tation (a — 1), and r — 20 gravitational radii, one 
can find (from Table 1) r = 1.8, which in Fig. 1 cor- 
responds to time interval T — t = 2itMc~ 1 t = 1.5 hr. 
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Twin peaks separation in sources with kHz QPOs caused by orbital 

motion 2 

Vladimir Karas 1 



The expected range of kilohertz quasi-periodic oscillations (kHz QPOs) is studied regarding low-mass X-ray binaries with 
twin peaks in their Fourier spectra. The twin peaks are interpreted as combination of epicyclic and precesion frequencies 
in Schwarzschild metric. If the pairs are caused by general relativistic orbital motion of gaseous clumps around a neutron 
star, then moderate size of the clumps and eccentricity of their orbits can account for the non-constant peaks separation 
reported in Scorpius X-l and 4U 1608-52. The neutron stars mass comes out in the range 1.6-1.9 solar masses. Variation 
in the orbital parameters of the clumps translates into an expected range in the plot of pair separation versus frequency, 
typically showing a decreasing slope. 
Keywords: X-rays: stars — stars: individual (Sco X-l, 4U 1608-52) — relativity 
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7. Introduction 

Accreting low-magnetic field neutron stars exhibit com- 
plex behaviour including X-ray variability on millisecond 
time-scales; see van der Klis (1998) for a review of mil- 
lisecond phenomena in low-mass X-ray binaries, and its 
recent update (van der Klis 1999). Cf. Mendez (1999) for 
a recent summary of RXTE observational results and for 
further references. Here we concentrate on the widely dis- 
cussed twin kHz peaks which so far have been reported in 
the Fourier power spectra of nearly twenty objects. The 
importance of this feature follows from the fact that the 
positions of the two peaks (and of the third, burst oscilla- 
tion frequency), their large coherence and variation with 
the count rate restrict viable models of the accreting body 
and of the geometry of the system. Also, it turns out that 
effects of general relativity cannot be ignored (Kluzniak 
1998; Stella & Vietri 1998; Zhang et al. 1998). 

The approximately constant separation of the twin 
QPO peaks indicates that some kind of beat-frequency 
relation (Miller, Lamb, & Psaltis 1998) can explain these 
pairs, although the simplest models have been challenged 
by the discovery of QPO sources in which the separation 
of the peaks decreases as their frequency increases: Sco 
X-l (van der Klis et al. 1997) and 4U 1608-52 (Mendez 
et al. 1998). Other sources are consistent with this trend 
too, e.g., 4U 1735-44 (Ford et al. 1998) and 4U 1728-34 
(Mendez &: van der Klis 1999). An inspiring alternative 
has recently been proposed by Stella & Vietri (1999) who 
introduce frequencies of the orbital motion around a com- 
pact body, namely orbital frequency, epicyclic frequency, 
and their difference. Here, this idea will be further stud- 
ied. In particular, the range of expected frequency varia- 
tion of the twin peaks will be estimated assuming that the 
feature arises from the motion of clumps with small but 
finite size, orbiting along eccentric orbits. The aim is to 
examine the main difference of this interpretation against 
simple constant-separation scenario by quantifying the ex- 
pected separation as a function of peaks frequency. The 
mass of the neutron star stands as a parameter of the pro- 
cedure in which predicted frequencies are fitted to data. 
It will be shown that the picture is indeed in agreement 
with the decreasing difference of the QPO peaks in Sco 
X-l and 4U 1608-52. Even more importantly, the model 
is testable in principle, and it can be rejected when pre- 
dicted frequencies of the peaks are found inconsistent with 
more accurate data in future. 

8. Frequencies relevant for the twin peaks 

Gravitational field is described in terms of Schwarz- 
schild metric (specified by the mass M of the body which 
is taken to be of the order of two solar masses). We will 
first write down expressions for the two fundamental fre- 
quencies driving the orbital motion of blobs, which are 



assumed to follow free trajectories (Chandrasekhar 1983). 
This is the simplest situation which of course will have to 
be abandoned in a more realistic model when the influence 
of radiation and of magnetic forces is taken into account, 
or when the gravitational field itself contains higher mul- 
tipoles. Trajectories in a spherical gravitational field of 
Schwarzschild metric are planar but, in general, they do 
not form a closed curve in space. In the case of bound 
stable orbits satisfying the condition 
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6^i-2^e>0, n,=Ma ^l-e 2 )" 



(1) 



one can associate eccentricity e and semimajor axis a with 
corresponding values of pericentre (r p ) and apocentre (r a ) 
distances, e = (r a —r p )/(r a +r p ), a — r p /(l — e). The orbit 
can be characterized by the pericenter shift 8<f) p during pe- 
riod P of time which elapses between subsequent passages 
through pericentre of the orbit. (Here, temporal quanti- 
ties are expressed with respect to a distant observer at 
rest and measured in geometrized units, while azimuthal 
angles are defined in the orbital plane in the standard 
manner). Both 5<f> p and P can be written in the form of 
integrals over relativistic true anomaly \, 



[l-2 M (3 + 7 7 )]- i/ "dx-27r 



(2) 



M (l + r,)- 2 [l-2 M (3 + r / )]- 1 / 2 

Jo 

x{l-2 f i(l + vT 1 dx, (3) 



where rj = ecosx, M = M\x 3/2 ^/(2^ - 1) - (2/ie) 2 . 
The frequency of radial oscillations (epicyclic frequency) is 
thus v r — P -1 , and precession frequency v p — 8<j> p /{2-kP). 
Combination of these two frequencies gives the corre- 



sponding azimuthal frequency, v$ 



Eqs. (2)- 



(3) can be written in terms of elliptic integrals (Karas & 
Vokrouhlicky 1994). Frequencies are then given in a more 
elegant form, although direct numerical evaluation of in- 
tegrals (2)-(3) is equivalent for practical purposes. 

Several instability mechanisms have been identified 
which may set clumps of gaseous material onto an eccen- 
tric orbit down from the inner edge of an accretion disk 

3 These quantities are introduced in standard manner, sim- 
ilar to Stella & Vietri (1999). Notice, however, that our 
definition of v p differs slightly from their corresponding re- 
lation for Vper, because here it is expressed in terms of the 
integral over the whole range of 27r in true anomaly. This 
is a proper and natural definition although the difference 
is only moderate for low-eccentricity orbits, considered in 
these papers. For eccentric orbits around 6M, azimuthal 
and precession frequencies are comparable while frequency 
of radial oscillations is much less. 
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Fig. 1. — The best fits in the plot of the pair separation 
Av(u) for v = z/0 (dot-dashed) and v — v v (dashed), re- 
spectively. For the latter case, shaded regions indicate the 
expected variation Av as function of the upper peak fre- 
quency v. Notice that interpretation of peaks in terms of 
epicyclic and precession frequencies with reasonable vari- 
ation of free parameters covers the data and captures the 
trend of decreasing Av{y). Frequencies are in Hertz. See 
the text for further details. 



or from the sonic point. Miller et al. (1998) remark that 
there may even be hundreds of clumps accumulated near 
the sonic point which can subsequently be set on an inspi- 
ralling trajectory. There is an interesting possibility which 
concerns freely moving clumps: a steady eccentric ring 
may be formed from individual blobs, surviving a number 
of orbital periods and rotating with precession frequency 
v v (Bao, Hadrava, & 0stgaard 1994). This presumes that 
clumps of matter are set on the eccentric orbit (a, e given) 
near the apocentre, one after another. Radiation from il- 
luminated or self-radiating clumps is modulated at orbital 
and precession frequencies, and their differences also ap- 
pear in the harmonic content of the observed signal. We 
thus examined the possibility of interpreting the differ- 
ence between the two peaks as Av = i/ r , while the upper 
twin peak can be interpreted either as v — v$ (as in Stella 
& Vietri 1999), or alternatively as v = v v . As we show 
later on, both interpretations lead to acceptable agree- 
ment with data on the peak frequencies in Sco X-l and 
4U 1608-52, but the latter one corresponds to the masses 
of the neutron star which are lower by about 0.3Mq. 

The frequencies were evaluated for orbits satisfying 
condition (1). (Notice that pericenter of eccentric orbits 
can be below the marginally stable orbit of circular mo- 
tion, which is at r — 6M in Schwarzschild metric.) There 
are three free parameters of the model, M, a, and e, which 
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Fig. 2. — Standard x 2 test of the best-fit mass parameter 
M (units of solar mass). In this figure, curves (i) and 
(ii) are for Sco X-l; (iii) and (iv) are for 4U 1608-52. 
Higher values of the mass, (ii) and (iv) (M — I.92M0), 
correspond to v — v<j,, while lower values, (i) and (iii) 
(M = I.6IM0), correspond to v = v v for both sources. 



can be varied to determine the best fit to the observational 
relation Av{y). Given these parameters, we still allow for 
a small variation of the integration range in x> Eqs. (2)- 
(3), which represents uncertainty in the clumps size d of 
the order of 3% of a. In other words, the bounds of the 
interval (0, 2ii) for true anomaly were varied by 8\{d) cor- 
responding to finite-size of clumps. For this purpose one 
only needs to know the relation between a small variation 
in true anomaly and the corresponding variation in az- 

imuthal angle: 8\ = £ [l — k 2 cos 2 (x/2)l 6<f>, where 
k 2 — 4/ie^" 1 and £ = 1 — 6/1+2/ie. Near apocentre, jfRJf, 
the shape of the trajectory is very close to its osculating 
ellipse in this point. 

Figure 1 shows the range which can be covered by this 
model in plots of Av{v) for Sco X-l and 4U 1608-52. 
These sources are known to exhibit decreasing separation 
Av when v increases (data taken from van der Klis et al. 
1997, and Mendez et al. 1998). The best fit is plotted 
for the two interpretations of the upper-peak frequency: 



(i) 



v,f, (dot-dashed line) and (ii) 



v p (dashed), 



respectively. We used standard \ statistics to determine 
the best-fit values of the parameters. The case (i) corre- 
sponds to Stella & Vietri (1999), and it yields the neutron 
star mass equal to 1.92M©, incidentally for both sources, 
in good agreement with their fit (apart from small differ- 
ence in calculation of relevant frequencies, which trans- 
lates to a 2% decrease of the best-fit mass; cf. footnote 2). 
Q probability is 10 -5 , defined in terms of F function and 
evaluated with deviations as in Fig. 1 (Press et al. 1992). 
The fit can be slightly improved by varying e, as sug- 
gested by Stella & Vietri, and one could also use a robust 
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fit. The case (ii) results in a graph of Av{v) which also 
decreases monotonically, but the slope is flatter. We will 
concentrate on this case since the former one has already 
been examined. The best-fit mass is now significantly less, 
I.6IM0 (Q — 0.4). The dark shaded region in Fig. la cor- 
responds to variation in 5\ around the best fit. The shad- 
ing covers a part of the (Av, u p ) plane where the best fit is 
moved by such variation. Although all parameters are in- 
terconnected in the equations for frequencies, one can say 
that M mainly determines the absolute value of upper- 
peak frequency v, while a and e affect the sloping stripe 
of Aviv). Typically, in our calculations a ranges between 
6M-10M, while r p can plunge down to AM, determining 
the corresponding eccentricity e < 0.4. Keeping a fixed, 
the slope increases as e decreases. An almost identical 
area of the graph can be covered by forcing the mass to 
a slightly higher value, M = 1.70Mq (indicated by some- 
what lighter shading in Fig. 1 for Sco X-l). In the latter 
case, however, 8\ has to be varied in a range correspond- 
ing to somewhat larger size of the clumps, from 3 to 6 per 
cent of a. For 4U 1608-52, the larger scatter in the data 
translates to a requirement on clumps size up to « 8% of 
a. 

Let us remark that non-zero size of the spots is in- 
evitable if they are to modulate visibly the observed radia- 
tion from the source for long periods of time. On the other 
hand, shearing motion destroys large spots and sets an up- 
per limit to their sizes. If individual elements of the spots 
followed free Keplerian trajectories on their neighbouring 
orbits, the above-mentioned 3% diameter (in units of o) 
would determine maximum spot size, which is still con- 
sistent with small observed widths of the peaks. 4 Indeed, 
the peaks' width increases with diameter of the sheared 
spots, and it should stay in agreement with the value of 
d, which has been determined here from scatter of data 
around the best fit in Fig. 1. Notice, however, that several 
physical mechanisms have been proposed where the spots 
in accretion disks could survive many orbital periods, in 
spite of their finite sizes. In this respect, formation of 
vortices was suggested by Abramowicz et al. (1992), and 
examined by Adams & Watkins (1995) and Bracco et al. 
(1999). It is speculated that these vortices durate over suf- 
ficient number of local orbital periods of the vortex center, 
giving rise to distinct features in observed power spectrum 
of the source. However, the origin of vortices is an open 
question so far. We thus remain on a phenomenological 
level of bright spots. 

Figure 2 shows graphs of x 2 , which were obtained by 
fitting a and e, while M was kept a fixed value within 



4 We thank the referee for this point. A particular example 
of a diminishing signal from spots which decay by shearing 
motion in the course of several revolutions was shown by 
Karas, Vokrouhlicky, & Polnarev (1992). 



the range 1.5-2.1 solar masses and 5\ = 0. Position of 
the minimum of x 2 (M) suggests the most probable mass. 
There are two different values of M clearly distinguished 
from each other, which correspond to the two cases of 
fitting described above, (i) M « 1.9M and (ii) M « 
I.6M0, although one must be cautious by interpreting 
such a graph when accuracy of data and the number of 
data points are limited. The fit can be further improved 
by factor ~ 2 by varying S\ as described above. 

Other physical processes, which have been neglected 
in the present text, introduce corrections to the expected 
range of frequencies, even if the main interpretation of 
relativistic orbital motion holds. Namely, one wants to 
estimate corrections which are caused by difference in the 
gravitational field of a rotating body in comparison with 
non-rotating Schwarzschild metric adopted here. It fol- 
lows from the discussion of Morsink & Stella (1999) that 
no realistic equation of state for the neutron star can ac- 
count for a change in kHz QPO frequencies more than 
10%. An analogous conclusion holds for motion in Kerr 
metric with maximum rotation. It is, however, impor- 
tant to note that inclined orbits around a rotating body 
are not planar. This is due to Lense-Thirring precession 
which may give rise to another fundamental frequency 
from a tilted accretion disk (Stella & Vietri 1998; Cui, 
Zhang, & Chen 1998). We refer to Merloni et al. (1999) 
for detailed evaluation of these fundamental frequencies 
of spherical orbits, and for corrections to Cui et al. (1998) 
paper. Even in the case of non-spherical perturbations 
of the orbits, which should occur in realistic situations, 
one can introduce a mean pericenter shift as an average 
taken over many orbits, or, equivalently, over a single or- 
bit of many individual clumps described by the proba- 
bility distribution of their orbital parameters. See Karas 
& Vokrouhlicky (1994) for such discussion of fundamen- 
tal frequencies of eccentric orbits in Kerr metric. Possible 
appearence of averaged frequencies in the power spectrum 
is illustrated there too. Notice that small non-zero eccen- 
tricity of the orbits of the clumps gives rise to variations 
of the expected peaks' frequencies about the mean value. 

9. Conclusion 

We have verified that observed frequencies of the twin 
kHz QPO peaks can be reproduced by assuming modu- 
lation of the observed signal by the two fundamental fre- 
quencies of orbital motion. An alternative estimate of the 
neutron star mass in Sco X-l is consistent with M ~ (1.6- 
1.7)M©, somewhat less than w 1.9Mq derived by Stella 
& Vietri (1999). The difference arises from interpreta- 
tion of the upper peak (y v vs. v$), and both possibilities 
still remain consistent with observational data. For orbits 
with negligible eccentricity we obtained identical results 
as those illustrated in their Fig. 1, but the present scheme 
obviously requires non-zero eccentricity to show up the 
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precession frequency in the signal. In both cases the ex- 
pected range of the pair's frequency corresponds to orbital 
motion around a < 2Mq compact object, and their sepa- 
ration decreases with increasing frequency. The decrease 
is faster for less eccentric orbits (within the relevant range 
of frequencies, i.e. 700 Hz < v < 1200 Hz for the upper 
peak). An appealing feature in the present highly sim- 
plified treatment is that the expected range, as indicated 
by the shading in Fig. 1, is very sensitive to the value of 
the central mass, and it moves away from the scale of our 
graphs when M is changed by ~ 0.3M & , independently 
of other free parameters. Although the orbital parame- 
ters e and 5\ are not strictly determined in the present 
treatment, it is worth mentioning that one cannot fiddle 
their values to obtain an increasing slope of Av(v); this 
might be possible only if v <C 700 Hz and/or M <J 1.5Mq. 
This argument can be used to exclude the present inter- 
pretation in other sources. For example, in the case of 4U 
1728-34 (Mendez & van der Klis 1999) the available data 
do not allow a consistent fit by the above described pro- 
cedure, although the peaks separation is again decreasing 
function of frequency. 
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Quasi-Periodic Features due to Clumps Orbiting around a Black Hole" 
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The expected Fourier power-spectrum density is examined with regard to the shape and quality factor Q of features which 
are predicted to appear in a signal from clumps orbiting near a marginally stable orbit of a compact body. It is shown that 
clumps distributed randomly in a narrow range of radii produce quasi-periodic oscillations with Q « 10 2 , similar to those 
observed in many low-mass X-ray binaries. 
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10. Introduction 

Different types of astronomical objects exhibit periodic 
modulation of their signal, which contains information 
about the physical properties of the source. Accreting 
low-mass X-ray binaries (LMXRBs) with quasi-periodic 
oscillations (QPOs) represent a particular example which 
has attracted much interest (for recent review see van der 
Klis 1999). While the origin of the QPOs is not fully 
understood, virtually all viable models introduce inho- 
mogeneities (clumps of gas) in an accretion disk orbiting 
near to a compact neutron star or a black hole. The fea- 
tures which are clearly visible in the Fourier power spectra 
can be attributed to Doppler and lensing effects modulat- 
ing the observed signal which comes from those rapidly 
moving clumps. Oscillations also arise due to some kind 
of beat frequency with neutron star spin (Miller et al. 
1998). The QPO frequencies can be rather high, up to 
~ 1300 Hz, which is comparable with dynamical time- 
scale at the marginally stable orbit of a stellar-mass ob- 
ject, Rms, suggesting distances of only several gravita- 
tional radii from the central body. Similar oscillations 
have been reported in X-rays from some active galactic nu- 
clei (AGN; Papadakis, Lawrence 1995). Here, the central 
masses are higher by six orders of magnitude and more, so 
that dynamical time is a few hours. Observational data 
are somewhat scarce in the case of AGN, compared to 
Galactic sources, because the time-scales are longer and 
the corresponding features in the power spectrum are less 
sharp than in the case of LMXRBs. Still, it appears that 
orbiting inhomogeneities contribute to the observed X- 
ray variability of AGN although their presence does not 
explain the entire variability spectrum on a purely kine- 
matic basis. The case of AGN is an issue of much debate, 
and other (magnetohydrodynamic) processes are clearly 
important (Abramowicz et al. 1991; Karas 1997). 

We examined the expected variability power spectrum 
from the distribution of clumps orbiting around a compact 
body. In this note we employ a phenomenological descrip- 
tion of the clumps, which has been adopted, in some form, 
by present models of sources with QPOs. These inho- 
mogeneities are visible in a narrow zone of a few gravi- 
tational radii (where luminosity of the accretion disk is 
maximum) with the lifetime greatly exceeding the corre- 
sponding orbital periods. The above-mentioned assump- 
tion concerning the strictly localized radial distribution of 
the clumps is one of the distinctions of the present work 
from analogous studies in which the short-term featureless 
X-ray variability of AGN has been explored in terms of a 
clumpy accretion disk. Bao & Ostgaard (1995) assumed 
a broad range of radii where the clumps exist in the disk 
and modulate observed signal (of the order of 10 2 gravita- 
tional radii), resulting in a power-law form of the Fourier 
power-spectrum density; in these circumstances individ- 
ual peaks are smeared away. The fundamental frequency 




Fig. 1. In this model, clumps are distributed 
in a zone of several gravitational radii 
with their observed luminosity influenced 
by Doppler and lensing effects. 



of orbital motion would stand out if all of the clumps were 
aligned along a ring with a strictly determined radius. Be- 
cause this will not be the case in more realistic models, 
it appears relevant to ask how the width of the feature 
in the power spectrum is modified when the clumps are 
spread over a small, but non-zero, range of radius. There 
are at least two reasons why a localized distribution of 
clumps should be taken into account. First, several insta- 
bility mechanisms have been proposed which could oper- 
ate within distances of several R g from the center of rela- 
tivistic accretion discs and result in trapped perturbations 
there (Kato, Fukue 1980; Markovic, Lamb 1998). Second, 
a non-axisymmetric structure of the disk surface could be 
best visible from the innermost parts of the disk where its 
radiation and irradiation are much stronger than farther 
from the center (e.g., Matt et al. 1991; Pariev, Bromley 
1998). Peaks in the Fourier spectrum corresponding to 
the localized distribution of the clumps can be charac- 
terized by the quality factor Q — v/Av (peak's centroid 
frequency divided by its full width at half maximum), and 
by harmonic content of the signal. 

We considered different distributions of the clumps, as 
characterized by their number density and local emissivity. 
Here we report the main conclusion of this study, namely, 
we show that maximum expected values of parameter Q 
are of the order of one hundred and they can be as large as 
3 x 10 2 . We will also show typical examples of the expected 
power spectrum with its harmonic content which is rather 
strong at medium and large (i.e. edge-on) inclination. 

11. Q-Factor and Harmonic Content of the 
Power Spectrum 

We assumed modulation of the observed signal by self- 
radiating or irradiated clumps orbiting on circular orbits 
in the Schwarzschild metric. The superposition of inde- 
pendent contributions from a large number (N S> 1) of 
such clumps and the noise component are the only sources 
of variability in this toy model. A similar situation has 
been discussed in various contexts. For example, in the 
beat-frequency model of QPOs one considers modulation 
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Fig. 2. Logarithmic plot of the expected 
power-spectrum density. The two pan- 
els correspond to: (a) pure superposition 
of light curves of mutually independent 
sources orbiting just above the marginally 
stable orbit; (b) an additional contribu- 
tion of the high-frequency noise. Ver- 
tical bars indicate Keplerian frequencies 
at the edges of the clumps radial dis- 
tribution, i.e. the range of fundamental 
frequencies (^k(-Rmin), ^k(-Rmax))- Here, 
the frequency (horizontal axis) is in ge- 
ometrized units, while the power spec- 
trum (vertical axis) has an arbitrary scale. 
See the text for details. 



of the signal by clumps near the sonic point (Miller et al. 
1998). However, the assumption about circular orbits is 
clearly a simplification which calls for further discussion. 
Stella and Vietri (1999) suggest that eccentric orbits of 
the clumps could explain non-constant separation of twin 
kHz peaks detected, e.g., in Sco X-l and 4U 1608—52. 
Each clump represents a spot revolving with Keplerian 
frequency vy_ — M 1 ' 2 r~' i ' 2 /(2iv) (geometrized units) . (For 
illustration: Keplerian frequency in Schwarzschild metric 
with M = 2Mq ranges down from its maximum value of 
1.1 kHz which is acquired at i? ms .) All of these spots re- 
side in a plane which coincides with an equatorial disk. 
The propagation of light from a source near a black hole 
to a distant observer has been treated by many authors 
within different contexts; we used one of the codes devel- 
oped in the past to calculate the observed profiles of fre- 
quency integrated light curves within a geometrical-optics 
approximation (Karas et al. 1992; see Kato et al. 1998 for 
a recent text-book exposition of the subject and further 
references). The orbital radius i? s and the observer's in- 
clination 9q are two free parameters, which determine the 
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Fig. 3. As in figure 2, but for the radial dis- 
tribution of sources with two maxima at 
slightly different radii in the disk [maxi- 
mum emitted flux is superposition of two 
Gaussian curves, each given by equation 
(2)]. Two neighbouring peaks are de- 
termined by the Keplerian frequencies at 
radii corresponding to the local maxima 
in the number density of the clumps dis- 
tribution. 



profile produced by a single source. It turns out that 
the size d of the spot and its detailed shape are less im- 
portant provided d^R^; this condition is assumed here- 
after. Given the two parameters, numerically calculated 
light curves were approximated by simple formula which 
captures the shape of the profile (the Doppler and the 
lensing enhancements included). This approach enables 
us to carry out fast superposition of hundreds of mutu- 
ally independent sources spread over a certain range of 7? s 
(Karas 1997). To do this, one also adopts that lifetime of 
the sources is much longer than their orbital period; one 
must prescribe the radial dependence of the local emis- 
sivity (e.g., a power-law). The resulting light curve in the 
time domain is the sum F(0q) — X^ =i F n (R s ,6o), and the 
power spectrum, P{v), is obtained by the Fourier trans- 
form F— >_F in the standard manner (Born, Wolf 1964): 



IP = linu 



t~ L \F\ 



The clumps are all placed slightly above the marginally 



stable orbit, R n 



3R g (where gravitational radius in 



terms of the central mass is R g = 2.95x 10 M/Mq). Their 
distribution is characterized by the range of radii AR be- 
tween R min and i? max (7? m in>-Rm S ; see figure 1), and by 
the corresponding radial dependence of the maximum flux 
from individual clumps, 



F n , 



oc R, 



(1) 
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where ai = const. We assumed locally isotropic emis- 
sivity of all sources in the frame corotating with the disk 
matter. The clumps are generated with random orbital 
phases and discarded after several orbits. The fluctuating 
signal is mean-subtracted and submitted to the Fourier 
analysis. 

We also explored various modifications of the clumps 
distribution; particularly, the power-law (1) (with sharp 
edges at i? m in and i? m ax) can be substituted by a smooth 
Gaussian-type profile, 



F n 



oc exp [— 011 [Rs — Re 



(2) 



where R ca n — Rmin + |Ai?, Q2 = const > 0. Such a flux 
distribution leads to a more familiar shape of peaks in the 
power spectrum, but, obviously, the choice of a particular 
relation requires one to specify the underlying physical 
mechanism producing and maintaining the clumps (which 
is beyond the scope of the present paper). It has been 
proposed and discussed by several authors that a gap in 
the disk can be opened where gas is swept out. This is 
caused by resonances due to the commensurability of the 
orbital period of gas with a distant disturbing secondary, 
or by the direct presence of a dense body (a 'planet') in 
the disk, which provokes the formation of inhomogeneities 
on the boundary of the gap. The resulting distribution of 
spots can then be approximated by superposing several 
profiles (2) with different i? cen and <xi 's. 

Figure 2a shows the power-spectrum density predicted 
by our simulation. The graph exhibits a typical peak 
around the fundamental orbital frequency of the clumps, 
v ~ (3-4) x 10 _3 cm _1 , and its harmonics. This plot corre- 
sponds to eq. (2) with N = 250, R min = 11M, AR = 3M, 
#o = 75°, and Q2 = 6. Conversion to physical units is 
achieved by multiplying the frequencies in geometrized 
units by the factor 2 x I0 5 (M/M ) _1 . In figure 2b, in 
order to account for high-frequency noise, we also super- 
posed simulated data with a noise component character- 
ized by zero mean and a standard deviation of a = 5%. 
The noise is naturally enhanced at higher frequencies, but 
does not greatly influence the values of Q, which refer to 
the peaks in the Fourier spectrum at lower frequency. 

Figure 3 illustrates a more complicated situation when 
the radiation flux emitted by the clumps in the disk has 
two maxima at different radii (superposition of two pro- 
files (2) with 9u = 25° and different values of i?cen,i = 
15M, i? C on,2 = 17 M; other parameters as in figure 2). Mo- 
tivation for such double-peaked radial distribution arises 
from possible formation of a gap in the disk. The funda- 
mental peak in the resulting power spectrum is split into 
two neighbouring components, and the signal has a strong 
harmonic content even at moderate inclination. 

Figure 4 shows the resulting quality factors for different 
distributions of the clumps. Here, we give the maximum 
expected Q corresponding to equation (1) with no noise 
component present in the signal (the values are decreased 
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Fig. 4. Quality factor Q plotted as function of 
the dimensionless radius i? m i n /M (mini- 
mum radial distance of the clumps) . This 
plot illustrates how Q increases with i? m i n - 
The values of Q are less if the clumps are 
distributed over wider range of radii, as 
characterized by Ai?. The three lines cor- 
respond to AR = \M (solid); AR = 1M 
(dashed); Ai? = 3M (dotted). Q is only 
weakly sensitive to other parameters (see 
the text). 



by noise, which increases the width of the peaks). Other 
parameters of the distribution were N = 250, ai = 1. 
Here, the values of Q were averaged over different incli- 
nations: 6o = 25°, 50°, and 75°. It turns out that Q 
depends only weakly on 0q. This conclusion can be under- 
stood by realizing that the redshift factor depends (up to 
relativistic corrections; Cunningham 1975) on cos^o and 
determines both the centroid energy of radiation from an 
orbiting source and the full width of the observed feature. 
Consequently, the #o-dependence of Q is not particularly 
strong provided one accepts the restriction 9q < 75° and 
r^,Rms, as we do in this graph. However, let us note that 
the harmonic content of the signal does depend on the in- 
clination. This is because the light curve of a source orbit- 
ing at small inclination resembles a sine profile, while it be- 
comes distorted by the lensing effect for large #o (^ 70°). 
In other words, the power in higher harmonics increases 
when the inclination increases and/or the mean lifetime of 
the spots decreases (to be specific, we assumed 5 orbital 
periods as typical lifetime in our examples). The power 
spectrum thus has potential to reveal orientation of the 
disk plane. On the other hand, an excessive amount of 
power in the overtones can be used to restrict the validity 
of the model. The above-given results are not very sensi- 
tive to the assumed local size of the clumps, provided that 
it is much less than R s , and to the life-time of the spots 
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provided it is much longer than corresponding orbital pe- 
riod. 

In our toy model, it appears quite natural that Q 
increases with the radius, as can be seen in figure 4 
(AR/Rmi n decreases when radius increases). One should, 
however, notice that Q > 10 2 also near i? ms . These results 
are not very sensitive to the values of the free parameters, 
and one can thus consider comparisons with observational 
data for particular objects (work in progress). 

12. Conclusions 

We presented a quantitative study of the expected 
variability power spectrum relevant for models of quasi- 
periodic oscillations. The orbital motion of the clumps 
and their distribution within a narrow range of radii are 
essential. The adopted phenomenological description of 
the source employs clumps orbiting in the Schwarzschild 
metric, and it can be considered to be the first approxima- 
tion to more realistic models of inhomogeneous disk-type 
accretion flows around black holes and neutron stars. Al- 
though incompleteness of the present treatment results in 
free parameters which are not strictly determined here, we 
checked that the main results (namely, the quality param- 
eter Q) are not strongly sensitive to reasonable variations 
of parameter values. 

The discussion of the power spectrum presented in this 
paper is a standard starting point of a Fourier analysis. 
It will be extended further by also examining the Fourier 
phases and their possible, coherence which should natu- 
rally appear due to the orbital motion of the clumps with 
the lifetimes substantially exceeding the dynamical time 
(Abramowicz et al. 1997; Karas 1997). It appears that the 
wavelet analysis (e.g., Starck et al. 1998) is better suited 
for this task than more usual approach which is offered by 
Fourier analysis. 
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